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Abstract:  A  new  model  is  derived  and  analyzed  for  tropical-extratropieal  inter¬ 
actions  involving  the  Madden-Julian  oscillation  (MJO).  The  model  combines  (i) 
the  tropical  dynamics  of  the  MJO  and  equatorial  baroclinic  waves  and  (ii)  the 
dynamics  of  barotropic  Rossby  waves  with  significant  extratropical  structure, 
and  the  combined  system  has  a  conserved  energy.  The  method  of  multiscale 
asymptotics  is  applied  to  systematically  derive  a  system  of  ordinary  differential 
equations  (ODEs)  for  three- wave  resonant  interactions.  Two  novel  features  are  (i) 
a  degenerate  auxiliary  problem  with  overdetermined  equations  due  to  a  compat¬ 
ibility  condition  (meridional  geostrophic  balance)  and  (ii)  cubic  self-interaction 
terms  that  are  not  typically  found  in  three- wave  resonance  ODEs.  Several  exam¬ 
ples  illustrate  applications  to  MJO  initiation  and  termination,  including  cases  of 
(i)  the  MJO,  equatorial  baroclinic  Rossby  waves,  and  barotropic  Rossby  waves 
interacting,  and  (ii)  the  MJO,  baroclinic  Kelvin  waves,  and  barotropic  Rossby 
waves  interacting.  Resonance  with  the  Kelvin  wave  is  not  possible  here  if  only 
dry  variables  are  considered,  but  it  occurs  in  the  moist  model  here  through 
interactions  with  water  vapor  and  convective  activity. 
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1  Introduction 

The  Madden-Julian  Oscillation  (MJO)  is  the  dominant  component  of  intrasea- 
sonal  («30-60  days)  variability  in  the  tropics  [10,  11,  12],  It  is  an  equatorial  wave 
envelope  of  complex  multi-scale  convective  processes,  coupled  with  planetary- 
scale  («10, 000-40, 000  km)  circulation  anomalies.  Individual  MJO  events  prop¬ 
agate  eastward  at  a  speed  of  roughly  5  m/s,  and  their  convective  signal  is  most 
prominent  over  the  Indian  and  western  Pacific  Oceans  [39].  In  addition  to  its 
significance  to  its  own  right,  the  MJO  also  significantly  affects  many  other  com¬ 
ponents  of  the  atmosphere-ocean-earth  system,  such  as  monsoon  development, 
intraseasonal  predictability  in  mid-latitude,  and  the  development  of  the  El  Nino 
southern  oscillation  (ENSO)  [8,  39]. 

In  addition  to  its  strong  tropical  signal,  the  MJO  interacts  with  the  global 
flow  on  the  intraseaonsal  timescales.  Telecormection  patterns  between  the  global 
extratropics  and  the  MJO  were  described  in  an  early  observational  analysis 
by  Weickmann  (1983)  [36]  and  Weikmann  et  al.  (1985)  [38].  Their  results 
demonstrate  coherent  fluctuations  between  extratropical  flow  and  eastward- 
propagating  outgoing  longwave  radiation  (OLR)  anomalies  in  the  tropics.  In 
the  study  by  Matthews  and  Kiladis  (1999)  [21],  they  illustrate  the  interplay 
between  high-frequency  transient  extratropical  waves  and  the  MJO.  More  re¬ 
cently,  Weickmann  and  Berry  (2009)  [37]  demonstrate  that  convection  in  the 
MJO  frequently  evolves  together  with  a  portion  of  the  activity  in  a  global  wind 
oscillation. 

The  interactions  between  extratropical  waves  and  tropical  convection  have 
also  been  investigated  in  numerical  models.  To  view  the  extratropical  response 
to  convective  heating,  Jin  and  Hoskins  (1995)  [4]  forced  a  primitive  equation 
model  with  a  fixed  heat  source  in  the  tropics  in  the  presence  of  a  climatological 
background  flow  and  obtain  the  Rossby  wave  train  response  in  the  result.  To 
diagnose  the  more  specific  response  to  patterns  of  convection  more  like  those 
of  the  observed  MJO,  Matthews  et  al.  (2004)  [20]  forced  a  primitive  equation 
model  in  a  climatological  background  flow  with  patterns  of  observed  MJO.  The 
resulting  global  response  to  that  heating  is  similar  in  many  respects  to  the  ob¬ 
servational  analysis.  The  MJO  initiation  in  response  to  extratropical  waves  was 
illustrated  by  Ray  and  Zhang  (2010)  [27].  They  show  that  a  dry-channel  model 
of  the  tropical  atmosphere  developed  M JO-like  signals  in  tropical  wind  fields 
when  forced  by  reanalysis  fields  at  poleward  boundaries.  Many  other  interesting 
studies  on  tropical-extratropical  interactions  have  been  carried  out.  For  exam¬ 
ple,  see  Lin  et.  al.  (2009)  [9]  and  Frederiksen  and  Frederiksen  (1993)  [3]  ,  and 
the  review  by  Roundy  (2011)  [32], 
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How  can  one  model  the  two-way  interaction  between  MJO  and  extratropi¬ 
cal  waves  in  a  simplified,  integrative  way?  This  is  the  primary  question  of  the 
present  paper.  Section  2  introduces  a  planetary  scale  model  for  this  purpose.  The 
model  includes  (i)  barotropic  dynamics  that  span  the  tropics  and  extratropics, 
(ii)  equatorial  baroclinic  dynamics,  and  (iii)  the  interactive  effects  of  moisture 
and  convection.  More  specifically,  the  model  integrates  the  dry  barotropic- first 
baroclinic  interaction  that  has  been  studied  by  Majda  and  Biello  (2003)  [14]  and 
Khouider  and  Majda  (2005)  [5]  with  the  MJO  skeleton  model  first  developed  by 
Majda  and  Stechmann  (2009)  [18]  and  further  investigated  by  same  authors  [19]. 
The  term  “skeleton”  used  by  the  authors  refers  to  the  fundamental  features  of 
the  MJO  on  planetary/intraseasonal  scales.  In  these  previous  studies,  Majda 
and  Biello  (2003)  carry  out  a  multiscale  asymptotic  analysis  for  the  resonant  in¬ 
teractions  of  dry  baroclinic  Rossby  waves  and  barotropic  Rossby  waves;  and  the 
MJO  skeleton  model  [18]  has  captured  three  main  features  of  the  MJO  on  plan¬ 
etary/intraseasonal  scales:  (i)  slow  eastward  phase  speed  of  roughly  5m/s;  (ii) 
peculiar  dispersion  relation  with  dco/dk  «  0;  (iii)  horizontal  quadrupole  vortex 
structure. 

A  multiscale  asymptotic  analysis  is  presented  here  in  Section  3  and  later 
sections,  adopting  similar  strategies  as  in  Majda  and  Biello  (2003),  based  on 
long-wave  scales,  small  amplitude  assumption,  and  multiple  long  time  scales.  The 
long-wave  scaling  leads  to  meridional  geostrophic  balance  at  the  leading  order, 
which  appears  as  a  constraint  in  the  partial  differential  equation  (PDE)  system. 
This  constraint  complicates  the  auxiliary  problem  of  the  multiscale  analysis.  In 
brief,  the  auxiliary  problem  is  a  necessary  step  for  suppressing  secular  terms.  The 
complication  here,  which  arises  from  the  constraint,  is  a  degenerate  operator, 
which  we  rectify  through  a  change  of  basis  as  described  in  Section  4.  Section  5 
presents  the  dispersion  relation  of  the  leading  order  linear  operator  for  both 
the  baroclinic  and  barotropic  modes.  In  Section  6,  by  inspecting  the  dispersion 
curves  of  the  linear  operator,  three-wave  resonant  interactions  are  identified,  and 
an  ODE  system  for  wave  interactions  is  derived  using  multiscale  asymptotics.  A 
novel  feature  of  this  ODE  system  is  the  presence  of  cubic  self-interactions  terms, 
which  are  not  typically  found  in  three- wave  interaction  ODEs  (e.g.,  [1,  13,  23,  30, 
31]).  Here  the  cubic  self-interaction  arises  from  the  new  nonlinearity  in  the  MJO 
skeleton  model  involving  water  vapor  and  convective  activity.  Many  previous 
studies  have  been  presented  on  resonant  interactions  of  atmospheric  waves,  and 
they  usually  focused  on  “dry  dynamics”  (e.g.,  [2,  15,  24,  25,  26,  28,  29]).  See 
Khouider  et.  al.  (2013)  [6]  for  a  review. 

Readers  who  are  most  interested  in  physical  applications  can  skip  ahead 
to  Section  6,  where  the  reduced  ODE  model  is  presented.  The  results  with  the 
model  are  then  organized  as  follows.  First,  in  Section  7,  a  validation  study  is  pre- 
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sented  to  explore  the  time  scales  of  validity  of  the  asymptotic  model.  Second,  in 
Section  8,  numerical  simulations  are  presented  for  three-wave  interactions.  Two 
cases  of  three- wave  interactions  are  chosen  for  application  to  MJO  initiation  and 
termination  [27,  33]  and  tropical-extratropical  interaction.  Each  case  therefore 
involves  both  the  MJO  and  barotropic  Rossby  waves,  and  the  third  wave  is  ei¬ 
ther  a  dry  baroclinic  Rossby  wave  or  a  dry  Kelvin  wave.  Interactions  involving 
only  the  Kelvin  wave  and  the  barotropic  Rossby  waves  are  not  possible,  but  by 
including  the  MJO  mode  from  the  model  with  water  vapor  and  convection,  the 
Kelvin  waves  are  engaged  in  the  resonance  triad. 

The  details  of  the  multiscale  asymptotic  analysis  are  presented  in  the  ap¬ 
pendices.  In  Appendix  A  and  B,  a  meridional  truncated  system  is  formulated. 
Appendix  C  provides  the  explicit  formulation  of  the  auxiliary  problems.  Finally, 
the  details  of  the  derivation  of  the  reduced  ODE  system  are  given  in  Appendix  D. 


2  The  two-layer  equatorial  /3-plane  equations 


The  nondimensional  two-layer  equatorial  /3-plane  equations  for  the  barotropic 
and  baroclinic  MJO  skeleton  model  are  given  by 


^  +  v  •  Vv  +  yv1-  +  X7p  =  — 7/V  ■  (v  ®  v), 
V  ■  v  =  0, 

dv  _  |  _ 

—  +  v  ■  Vv  -  V0  +  wv  =  -v  ■  Vv, 
at 

f)f)  - 

—  +  v  ■  V<9  -  V  ■  v  =  52(Ha  -  Se ), 

'q+Q\7  -v=  -52(Ha-Sq), 


(la) 

(lb) 

(lc) 


(Id) 


(le) 


(If) 


These  equations  combine  the  MJO  skeleton  model  [18]  and  nonlinear  inter¬ 
actions  between  the  baroclinic  and  barotropic  modes  [14].  The  equations  have 
been  nondimensionalized  using  the  scales  listed  in  Table  1.  Here  v  =  (u,v)  and 
p  are  barotropic  velocity  and  pressure;  v  =  (u,  v)  and  6  are  baroclinic  veloc¬ 
ity  and  potential  temperature;  and  q  is  water  vapor  (sometimes  referred  to  as 
“moisture”).  The  tropical  convective  activity  envelope  is  denoted  by  S2a,  where 
5  is  a  small  parameter  that  modulates  the  scales  of  tropical  convection  envelope. 
Likewise,  52  is  also  applied  to  quantities  SB  and  Sq,  radiative  cooling  and  the 
moisture  source.  In  this  paper,  52 Se  and  52Sq  are  set  to  be  constants  for  energy 
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Par. 

Derivation 

Dim.  val. 

Description 

P 

2.3  x  10’11  m_1s_1 

Variation  of  Coriolis  parameter  with  latitude 

e0 

300  K 

Potential  temperature  at  surface 

9 

9.8  m  s~2 

Gravitational  acceleration 

H 

16  km 

Tropopause  height 

N2 

(g/0o)d8/dz 

o 

1 

Wl 

1 

to 

Buoyancy  frequency  squared 

c 

NH/tt 

50  m  s  1 

Velocity  scale 

xe 

V°/P 

1500  km 

Equatorial  length  scale 

T 

XJc 

8  hrs 

Equatorial  time  scale 

HN2e0/(jrg ) 

15  K 

Potential  temperature  scale 

H/tt 

5  km 

Vertical  length  scale 

H/(nT) 

0.2  m  s”1 

Vertical  velocity  scale 

c2 

2500  m2  s'2 

Pressure  scale 

Table  1.  Constants  and  reference  scales  for  nondimensionalization. 

Par. 

Non-dim.  val. 

Dim.  val. 

Description 

r 

1 

~0.5  /day/(g/kg) 

Convective  growth/decay  rate 

H 

0.23 

~10  K/day 

Parameter  to  rescale  a 

52a 

<52 

Convective  activity  envelope  at  RCE  state 

Q 

0.9 

Non-dim.  background  vertical  moisture  gradient 

82Se 

82H 

Radiative  cooling  rate 

82Sq 

52H 

Moisture  source 

Table  2.  Parameters  of  the  MJO  skeleton  model,  and  relation  to  small  parameter  8. 


conservation,  although  they  usually  have  both  spatial  and  temporal  variance  in 
reality.  Together  with  T,  H  and  Q ,  the  coefficients  are  described  in  table  2. 

In  equations  (1),  vx  =  (—v,u)  is  the  ,5-plane  approximation  of  tropical 
Coriolis  force,  and  x  and  y  denote  the  zonal  and  meridional  coordinates.  Without 
the  moisture  q  and  convection  envelope  a,  system  (1)  is  the  two-vertical-nrode 
Galerkin  truncation  for  the  Boussinesq  equations  with  rigid-lid  boundary  on  the 
top  and  bottom  of  the  domain  (see  e.g.,  [22,  16,  14,  5])  Without  the  barotropic 
wind  v,  the  system  is  the  MJO  skeleton  model  on  the  first  baroclinic  mode 
before  passing  to  the  long  wave  limit  (see  e.g.,  [18,  19]). 

Note  that  for  the  primitive  equations  (1),  a  total  energy  is  conserved,  and  it 
is  composed  of  four  parts,  dry  barotropic  energy  £?,  dry  baroclinic  energy  £q, 
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moisture  energy  Em ,  and  convective  energy  Ea- 

y  x 

£r(t)  =  \  J  J  lv|2dxdy  (2a) 

—Y  0 
Y  X 

Ec{t)  =  \  J  J  lv|2  +  02dxdy  (2b) 

-Y  o 

y  x 

+  <2c> 

-y  o 
y  x 

£a(*)  =  y  J  J  -jE  [Ha-  S6  log  (a)]  dxdy  (2d) 

-Y  0 

This  energy  conservation  forms  the  design  principle  as  appeared  in  [18,  14]. 
Note  that  in  (If),  the  conserved  quantity  £  still  holds  if  a  is  also  advected  by 
the  barotropic  wind,  i.e.,  with  an  additional  term  u  ■  Va  in  (If).  But  we  do  not 
include  this  term  because  later,  when  meridional  truncation  is  applied  to  the 
system,  the  energy  conservation  will  not  hold  with  this  additional  term. 

Using  the  streamfunction  tp  for  barotropic  mode,  which  satisfies  (u,v)  = 
(—ipyjipx),  the  barotropic  equation  can  also  be  written  as 

-^A^  +  ipx  +  ■  [-(vu)j  +  (vd)x]  =  0,  (3) 

where 

D  d  _  d  _  d 

-  = - b  U - b  V - 

Dt  dt  dx  dy 
represents  advection  by  the  barotropic  wind. 


2.1  The  zonal-long  wave  scaled  model 

To  consider  the  planetary  scale  of  the  coupled  equations  in  (3)  and  (lc)-(lf), 
zonal  variations  are  assumed  to  depend  on  a  longer  scale,  as  are  temporal  vari¬ 
ations.  The  long  zonal  and  long  temporal  coordinates  are  defined  as 

x'  =  5x,  t'  =  dt.  (4) 

Correspondingly,  the  meridional  velocity  is  also  scaled  so  that 


v  =  5v' . 


(5) 
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The  equations  in  (3)  and  (lc)-(lf)  become  the  long-wave-scaled  system: 


+  Vv  -  *  V'  ■  [(v'u)y]  +  S2  < 

J +  \v'  ■  [(v'nOx'] 

|=0,  (6a) 

■^7 u  -  6X  -  yu  -  V  ■  V'ipv  =  0, 

(6b) 

-6y+yu  +  S2  +  v'  ■  Vipx>^J 

=  0 

(6c) 

JLe  -  v  ■  v  -  5{Ha  -  Sd)  =  0, 

(6d) 

J^q  +  QV  -v'  +6(Ha-Sq)=  0, 

(6e) 

3?l<a 

p 

1 

n 

p 

II 

o 

(6f) 

Here  the  primes  represent  the  long-wave  scaled  coordinates  and  variables: 


D 

dv 


d  d  ,  d  ,  ,  d  d 

dt'  U  dx'  V  dy'  ’  dx'  ’  dy 


and  v'  =  (u,v'). 


System  (6)  is  in  the  same  form  as  in  [14]  when  the  moisture  and  convection  are 
neglected. 

Furthermore,  the  convective  activity  a  will  be  written  as  an  anomaly  from 
the  state  of  radiative-convective  equilibrium  (RCE): 

Se  =  Sq  =  Ha, 


where  these  constants  take  the  values  given  in  table  2.  When  a  =  a,  the  external 
radiative  cooling/moisture  source  is  in  balance  with  the  source/sink  from  the 
convection,  and  the  system  achieves  RCE.  By  writing  a  as  an  anomaly  with 
respect  to  the  RCE  value  a,  the  forcing  terms  in  (6d)  and  (6e)  can  be  written 
as 

Ha  -  Se  =  Ha  -  Sq  =  H(a  -  a)  =  Ha',  (7) 

and  the  equation  for  the  convection  envelope  (6f)  can  be  written  as 

■^-a'  -Taq  =  Ta'q.  (8) 


2.2  Long  time  scales  for  tropical-extratropical  interactions 

The  small  parameter  5  is  also  used  for  introducing  two  longer  time  scales: 

Ti  =  dt',  T2  =  52t'. 

Their  units  match  the  units  of  intraseasonal  timescales  as  appeared  in  [14,  17] 
In  the  next  section,  the  system  is  expanded  by  matching  the  orders  of  5. 
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3  Asymptotic  expansions  for  the  interaction  of 
barotropic  and  equatorial  baroclinic  waves 

In  this  section,  asymptotic  expansions  are  carried  out  twice:  first  for  the  2-D 
system  as  functions  of  x  and  y,  and  second  for  a  system  using  a  truncated 
basis  for  variations  in  the  y  direction.  The  first  expansion  (2-D)  provides  the 
clearest  presentation,  but  the  2-D  linear  operator  does  not  have  eigenvalues  and 
eigenfunctions  that  are  easily  accessible.  For  this  reason,  the  second,  truncated 
system  is  introduced  and  provides  a  linear  operator  with  known  eigenvalues  and 
eigenfunctions. 


3.1  The  2-D  equations 

Assume  that  the  solutions  have  an  asymptotic  structure  with  ansatz: 

(ip,u,v',6,q)  =  52(ipx,ux,vx,9x,qx)  +  53(ip2,  u2,  v2, 02, 42) 

+  d4{tl)3,u3,V3,03,q3)  +  0(S5)  (9a) 

and  a 1  =  Sax  +  52a,2  +  53ci3  +  0(5 4),  (9b) 

This  small  amplitude  assumption  is  consistent  with  a  small  Froude  number 
assumption  as  in  [14].  All  variables  are  assumed  to  depend  on  three  time  scales: 
t'  and  T\  =  5t'  and  T2  =  52t' .  The  system  is  expanded  over  three  orders  of 
magnitude. 

The  first  order  system  is 


4>X  y’y’t’  +  Iplx1  =  0,  (10a) 

wit'  -  Qix1  -  y'v  1  =  0,  (10b) 

-  9i y>  +  y’m  =  0,  (10c) 

6 w  ~  Mix'  -  viy>  -  Hax  =  0,  (10d) 

Qit'  +  Q(u  lx'  +  Vx  v')  +  Hax  =  0,  (10e) 

ait/  —  raq1!  =  0.  (10f) 


This  system  defines  the  leading-order  linear  operator,  including  a  combination 
of  barotropic  Rossby  waves  [14]  from  the  dynamics  of  ip  and  baroclinic  equatorial 
long-waves  including  the  MJO  [18]  from  the  other  variables. 
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The  second  order  system  is 

^2 y'y't'  1p2x'  —  V’ly'y'TY:  (lift) 

U2f  ~  fox'  ~  V’v 2  =  -MlT1;  (Hb) 

-  #2y'  +  y'u2  =  0,  (lie) 

02 1>  ~  U2x’  ~  V2 y>  ~  Ha2  =  -01Ti,  (Hd) 

<72t'  +  Q(u2x >  +  v2y' )  +  H a2  =  -qirx ,  (He) 

a2f  ~  ra?2  =  ~a1Tl  +Ta1qi.  (Ilf) 


Note  that  the  only  nonlinear  term  is  from  the  q-a  interaction:  Taigi. 

The  third  order  system  is 

i>3 y'y’t'  +  tp3x'  =  —Ipl y'y'T2  ~  ^ 2y'y'Tx 

—  Iplx'x't'  T  ^lx^ly'y'y1  ~  ^ly'^lx'y'y'  ~  ^  (^1  )x'y'  —  (t£l V\ 'jy'y'  j 

(12a) 

U3t'  —  93x'  -  y'v 3  =  -U\T2  —  U2Ti  —  [iplx'Ul y>  -  i>ly>Ulxi  -  U\lj)lx'v'  -  I’l'i/’l y'y']  , 


(12b) 

-  #3 y>  +  y'u3  =  -Wit',  (12c) 

03t'  u3xf  v3 y>  H&3  9it2  92T\  V’lx'^iy'  T  Vhy'  9\x'  j  (12d) 

<?3t'  +  Q{u 3x'  +  v3y')  +  Ha3  =  —qiT2  ~  q2Tx  -  'I’lx'qiy'  +  V’ly'Qix'i  (12e) 

a3f  ~  Faq3  =  -air2  —  d2Tx  +  r(aig2  +  a2qi)  —  ipix'diy1  +  ^ly'^ix'-  (12f ) 


Note  that  more  nonlinear  terms  arise,  including  the  baroclinic-barotropic  in¬ 
teractions  as  in  [14]. 

The  total  energy  is  approximately  conserved: 

^(3)  =  o(54).  (13) 

from  the  asymptotic  expansions  (10)-(12),  where  £ ^  is  the  energy  of  the  vari¬ 
ables  from  the  leading  3  orders  of  magnitude  from  the  ansatz  in  (9).  This  ap¬ 
proximated  energy  conservation  is  important  to  preserve  in  future  derivations. 

The  asymptotic  expansions  (10)-(12)  are  analogous  to  the  expansions  in  [14], 
where  the  system  did  not  have  moisture  q  nor  convection  envelope  a.  Nonethe¬ 
less,  there  are  essential  differences  described  in  table  3.  For  example,  in  [14], 
the  small  parameter  was  e  =  S2,  in  which  case  the  ansatz  in  (9)  includes  only 
the  52  and  S4  terms,  and  only  the  systems  in  (10)  and  (12)  arise.  On  the  other 
hand,  here  the  small  parameter  is  fi,  which  leads  to  the  additional  <53  terms  in 
the  ansatz  (9)  and  the  additional  system  in  (11). 
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MJO  skeleton  model 

Dry  wave  model  in  [14] 

Moisture  equations  involving  q  and  a  are  included 

Dry  dynamics  only 

Small  parameter  8 

Small  parameter  e  =  82 

5  =  •v/O.l,  e  =  0.1 

Two  long  time  scales:  7\  =  St'  and  Tb  =  52t' 

One  long  time  scale:  r  =  et' 

T\  =  0(3  day),  Tb  =  0(10  day) 

Linear  system  is  dispersive 

Linear  system  is  non-dispersive 

Table  3.  Comparisons  between  the  MJO  model  here  and  the  dry  wave  model  of  [14], 


Unfortunately,  the  2-D  linear  operator  defined  in  (10)  does  not  have  eigen¬ 
values  and  eigenfunctions  that  are  easily  accessible.  This  is  due  to  the  effects  of 
moisture  q  and  convective  activity  a,  since  on  the  other  hand  the  “dry”  dynamics 
without  q  and  a  do  have  well-known  eigenvalues  and  eigenfunctions  [14].  Since 
the  method  of  multiscale  asymptotics  relies  heavily  on  the  leading  order  linear 
operator,  further  progress  cannot  be  made  with  the  systems  ( 10) — (12). 

To  circumvent  this  issue,  a  truncated  basis  will  be  introduced  next  in  section 
3.2  for  variations  in  the  y  direction,  and  the  truncated  system  has  eigenvalues 
and  eigenvectors  that  have  previously  been  presented  [18].  The  truncated  system 
utilizes  Riemann  invariants  /  =  —(«  +  9)/y/ 2  and  r  =  (u  —  9)/y/2  [13]  in  place 
of  the  variables  u  and  9. 


3.2  Meridional  truncated  system 

To  simplify  the  system,  a  meridional  truncation  is  adopted  as  in  the  truncations 
in  MJO  skeleton  model  [18]  by  using  parabolic  cylinder  functions  for  baroclinic 
variables.  The  barotropic  variables  are  assumed  to  be  sinusoidal  as  in  [14].  Al¬ 
though  mostly  similar  to  [18]  and  [14],  some  modifications  are  needed  for  energy 
conservation,  as  described  further  below.  Explicitly,  the  meridional  trunca¬ 
tion  is 


ip  =  B(x ,  t )  sin (Ly) 

(14a) 

q)  =  (l{0\x,t),ri0\x,t),q('°\x,t) 

)  $o(2 /), 

(14b) 

+  (l(2)  (x,  t),r{2)  (x,  t ) ,  q{2)  (x,  t)' 

)  $2 (y), 

(14c) 

v'  =  v(1)4>i  (y), 

(14d) 

a'  =  a(0)(a;,f)$o(2/), 

(14e) 
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where  $m(y)  are  parabolic  cylinder  functions: 

®m(y)  =  (rralyV)  2  2 ^fe~^Hm(y),  (15) 


with  Hermite  polynomials  Hm(y)  defined  by 

2  rlmp~y2 

Hm(y)  =  (~l)mey  -1^r.  (16) 

The  parabolic  cylinder  functions  form  an  orthonormal  basis  on  the  ID  function 
space.  The  first  few  functions  are 


=  ^-^-2/72 


$0 (2/)  =  7r  4  e 


$i(y)  =  7r  !/2/2,  $2(2/)=' 


-l—(2y2-l)e~y2 

(17) 


The  parabolic  cylinder  functions  satisfy  the  following  identities: 

£+®m(y)  =  (2m)1/2$m_1(2/),  C-<bm(y)  =  -  [2(m+  1)]1/2  d>m+i(?/),  (18) 


which  help  to  simplify  the  expression,  where  the  operators  C±  are  defined  as 


C± 


d 

dy 


±y- 


In  the  truncation  (14a),  the  parameter 


L  = 


2t r 
2 Y 


is  the  meridional  wavenumber,  and  2Y  is  the  meridional  wavelength. 

Applying  this  meridional  truncation  to  the  2-D  model  in  (6)  leads  to  the 

truncated  system: 


L2YBt>  -  YBX ,  =Ct+52  {YBx,x,t,  +  VT) ,  (19a) 

MUt  +  C0U  =  Cc  +  52VGi  (19b) 

where  U  =  (l^°\  l^2\  A2\  fZ1)  ,q(°\q(2\a(  °) ) ,  is  the  spatial  linear 

operator,  and  J\f  =  diag(l,  1, 1, 1, 0, 1, 1, 1)  is  an  8  x  8  diagonal  matrix  with  the 
0  placed  in  the  element  of  the  diagonal.  At  the  right  hand  side,  Ct, 

Cc  and  T>q  are  bilinear  terms  at  different  orders  in  the  long-wave  scaled  system. 
These  bilinear  terms  are  from  advection  terms  in  the  2-D  system  (i.e. ,  v  ■  Vv, 
V  ■  (v  ®  v),  v  •  Vv,  v  ■  Vv)  and  nonlinear  q-a  interaction:  Tqa.  The  detailed 
expressions  are  given  in  (A.l)  in  Appendix  A. 

This  particular  meridional  truncation  (14)  is  used  because  it  maintains 
the  Z(°)_'i,(1)_r(2)  triplets  to  obtain  the  correct  dispersion  relation  for  baroclinic 
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Rossby  waves  [13].  Also,  this  meridional  truncation  imitates  the  energy  conser¬ 
vation  as  in  the  full  system  with  minimal  alteration  in  the  a ^  equation.  Other 
truncations  may  also  be  considered,  such  as  neglecting  l^>  and  q as  in  MJO 
skeleton  model  [18],  and  we  have  explored  several  other  options.  However,  other 
truncations  can  potentially  contribute  to  energy  imbalance  from  the  barotropic 
advection  acting  on  baroclinic  variables. 

Finally,  the  truncated  system  (19)  can  be  expanded  in  powers  of  5,  similar 
to  the  expansion  (10)-(12)  of  the  2-D  system  (6).  The  expansion  of  (19)  is  pre¬ 
sented  in  Appendix  B  in  (B.2)-(B.4).  The  next  step  in  the  method  of  multiscale 
asymptotics  is  to  solve  the  auxiliary  problem  that  arises  in  each  of  (B.2),  (B.3), 
and  (B.4). 


4  Auxiliary  problem:  a  degenerate  system  with 
constraint  equation 

In  this  section,  the  auxiliary  problem  is  introduced.  It  is  a  key  part  of  the  multi- 
scale  analysis  procedure  where  a  forced  linear  system  must  be  solved  (see  chapter 

5  of  [13]  for  other  examples).  In  the  linear  system  of  the  present  paper,  a  difficulty 
arises  because  of  a  constraint  equation;  specifically,  a  constraint  arises  from 
meridional  geostrophic  balance  where  does  not  appear  in  the  leading  order 
linear  operator  in  (10)  or  (A.l).  This  constraint  equation  presents  a  difficulty 
for  directly  computing  the  eigenmodes  of  the  linear  system.  To  overcome  this 
difficulty,  a  reformulation  of  the  auxiliary  problem  is  presented  here  to  enable 
direct  calculation  of  the  eigenmodes  and  hence  a  direct  solution  of  the  forced 
linear  system. 

The  linear  system  in  U  for  the  baroclinic  component  from  (19)  can  be  de¬ 
noted  as: 

Mtt  +  C&U  =  Fff.  (20) 

Notice  the  degenerate  diagonal  8x8  matrix  J\f  =  diag(l,  1, 1, 1, 0, 1, 1, 1),  which 
is  degenerate  due  to  the  0  in  the  uW  element  of  the  diagonal.  Here  F g  is 
the  forcing  term,  including  long-time  dependency  and  bilinear  terms.  Also  re¬ 
call  from  (19)  that  U  is  the  vector  of  baroclinic  variables,  defined  as  U  = 
(^(o),/(2),r(cl),r(2),i;(1),,g(0),gl(2),a(o))  pjue  -|-0  degenerate  matrix  M,  the 
eigenmodes  of  this  system  cannot  be  found  directly  using  the  standard  proce¬ 
dure. 
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In  order  to  circumvent  the  degenerate  matrix  J\f ,  the  auxiliary  problem  can 
be  reformulated  in  terms  of  new  variables: 

W=  (K,R,Q,A:v^,xJ{2\q(2)), 

where  K  =  R  =  \/2l^  +  2A2\  Q  =  q(°\  A  =  a^°\  x  =  ^  —  \/2 A2\  This 
reversible  linear  transformation  can  be  written  as  W  =  AU.  Then  the  system 
for  W  is 

ANA-'dtW  +  ACjjA-'W  =  AF0,  (21) 

which  can  be  further  simplified  to 

dtW  +  £^W  =  F^i  (22) 

where  W  =  (K,  R,Q,  A,l(2\qW),  and  is  the  spatial  linear  operator.  From 
W  to  W,  two  variables  and  \  are  eliminated  because  of  the  constraint 
equation.  For  details,  see  Appendix  C. 

Note  that  this  transformation  of  variables  from  U  to  W  is  used  in  the  for¬ 
mulation  of  the  MJO  skeleton  model  [18,  34],  and  it  has  some  similarities  to  the 
transformation  of  the  two-dimensional  incompressible  fluid  flow  equations  from 
velocity  variables  to  vorticity/streamfunction  variables  (see  chapter  9  of  [13]). 
For  the  purposes  here,  the  formulation  in  terms  of  W  =  (K,  R,  Q,  A,  q^) 
allows  easy  identification  of  the  linear  eigenmodes  of  the  linear  operator  £q>, 
since  the  (K,  R,Q,  A)  sub-system  of  (22)  is  the  linearized  MJO  skeleton  model 
that  has  been  previously  studied  [18] . 

Finally,  the  forcing  Fjj  in  (20)  must  be  transformed  into  F^  in  (22).  This 
transformation  is  needed  in  order  to  solve  the  auxiliary  problem  at  each  order 
of  magnitude  in  the  asymptotic  expansion.  See  Appendix  C  for  the  explicit 
formulas  for  F^  for  each  order  of  magnitude  of  the  asymptotic  expansion. 


5  Eigenmodes  for  the  linear  system 

In  this  section,  the  eigenmodes  are  described  for  the  baroclinic  system  and 
barotropic  system.  (Later,  in  Section  6,  these  linear  eigenmodes  will  be  used 
in  the  identification  of  three-wave  resonances.) 

First,  in  the  baroclinic  system  (22)  (with  F^  =  0),  the  components  of  vector 
W  =  (K,R,Q}A,l(2\q(2^)  can  be  separated  into  two  groups:  ( I\,R,Q,A )  and 
{l^2\q^).  The  KRQA  system  is  closed  (as  can  be  seen  in  its  explicit  formulation 
in  Appendix  C),  and  its  four  eigenmodes  have  been  described  previously  [18]: 
dry  Kelvin,  MJO,  moist  Rossby  and  dry  Rossby,  as  shown  in  figure  1.  In  brief, 
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Fig.  1.  Dispersion  relation  of  the  KRQA  system  and  barotropic  wind  with  different  merid¬ 
ional  wavelengths. 


the  dry  Kelvin  wave  is  a  fast  eastward-propagating  wave;  the  MJO  is  a  slow 
east-propagating  wave;  and  the  moist  and  dry  baroclinic  Rossby  waves  are  slow 
and  fast  westward-propagating  waves,  respectively.  In  addition  to  the  KRQA 
system,  l ^  satisfies  an  independent  (decoupled)  equation,  and  q ^  is  slaved  to 
( K ,  R,Q,A)  and  l ^  (see  the  explicit  formulation  in  Appendix  C).  Therefore  the 
eigenmodes  from  the  KRQA  system  of  [18]  are  not  affected  by  the  additional  two 
variables  l ^  and  q1'2'1 .  Notice  that  the  baroclinic  Rossby  and  Kelvin  waves  in 
this  model  are  dry  equatorial  waves,  not  convectively  coupled  equatorial  waves 
(CCEWs)  [7].  The  model  here  is  a  planetary-scale  model  that  does  not  explicitly 
resolve  CCEWs,  which  occur  mainly  on  smaller,  synoptic  scales.  Nevertheless, 
one  would  expect  similar  wave  interactions  to  hold  for  CCEWs,  due  to  the 
similar  structures  of  dry  and  convectively  coupled  waves,  if  a  model  were  used 
that  resolved  such  synoptic-scale  convectively  coupled  waves. 

Second,  for  the  barotropic  system  (19a),  the  dispersion  relation  is  written 
explicitly  as 

“T{.k)  =  ~.  (23) 

The  above  formula  implies  that  barotropic  Rossby  waves  travel  faster  for  smaller 
meridional  wavenumber  L  or,  equivalently,  longer  meridional  wavelength  Y .  In 
figure  1,  barotropic  dispersion  relation  are  shown  with  different  meridional  wave 
numbers. 
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6  The  reduced  asymptotic  models 


In  this  section,  the  following  reduced  asymptotic  model  is  derived  and  its  basic 
features  are  described: 


<9t2/?  +  id2/3  +  id^aXa^  =  0, 

dr2ai  +  id±ai2al  +  id$ai  +  id^jd*^  =  0, 

dr2ct2  +  id’ra22a2  +  idga2  +  idgfd*a\  =  0, 


(24a) 

(24b) 

(24c) 


<9t2/? 


where  the  coefficients  dj  are  all  real  numbers.  This  is  a  system  of  ODEs  for  three- 
wave  resonance,  where  /3  is  the  complex  amplitude  of  the  barotropic  Rossby 
waves,  and  Qi  and  «2  are  the  complex  amplitudes  of  two  baroclinic  waves.  The 
superscript  *  stands  for  complex  conjugate.  For  the  applications  of  importance 
here,  one  of  the  as  will  correspond  to  the  MJO,  and  this  ODE  system  describes 
its  interaction  with  the  other  waves,  with  an  aim  toward  MJO  initiation  and 
termination. 

In  brief,  the  terms  in  (24)  fall  into  three  groups  that  are  linked  to  different 
features  in  the  full  2D  model  (1):  the  cubic  terms  are  associated  with  the  nonlin¬ 
ear  q-a  interactions;  the  linear  terms  come  from  dispersive  terms;  the  quadratic 
terms  rise  from  the  nonlinear  baroclinic-barotropic  interactions.  Here,  in  partic¬ 
ular,  the  cubic  self-interaction  terms  are  a  novel  feature  that  are  not  typically 
found  in  three-wave  resonance  ODEs,  and  they  arise  here  from  the  effects  of 
water  vapor  q  and  convective  activity  a. 

In  the  remainder  of  this  section,  two  basic  properties  of  the  system  are 
described  -  energy  conservation  and  a  family  of  equilibrium  points  -  and  the 
derivation  is  sketched,  with  details  of  the  derivation  shown  in  Appendix  D. 

6.1  Energy  conservation 

The  ODE  system  (24)  satisfies  the  following  energy  conservation  principle: 


(25) 


This  conserved  quantity  is  a  consequence  that  all  coefficients  are  pure  imaginary, 
and  in  particular,  that 


d%  +  de  +  dg  —  0, 


(26) 


which  is  a  key  component  for  energy  conservation  in  three-wave  interaction 
equations  [13]. 
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6.2  Equilibium  points  in  ODE  system 

In  addition  to  the  trivial  equilibrium  point  of  (24)  where  /?  =  a\  =  a2 
family  of  nontrivial  equilibrium  points  also  exists: 

2  _  d\d3dy  +  d2d3deds 
d^dgdg  —  d^d^df 

p  i/-)  di  d3  +  d2d3dgd3 

d^dgdg  —  d^d/^d'j 

n  d3  *  jj, 

P  =  -q-ai  a2 
d  2 

Two  examples  for  the  nontrivial  equilibrium  points  are  considered  here 

|ai|w0.24,  \a2\  ~  0.20,  \ft\  &  0.0033  for  case  MRB,  (28a) 

|ai|f«2.4,  \a2\  ~  7.5,  \f}\  tv  0.34  for  case  MKB.  (28b) 

The  MRB  and  MKB  are  cases  for  different  resonance  triads.  For  MRB,  the  M 
stands  for  MJO,  the  R  stands  for  equatorial  Rossby  wave,  and  the  B  stands 
for  barotropic  Rossby  wave;  and  for  MKB,  the  M  and  B  are  the  same,  and  the 
K  stands  for  equatorial  Kelvin  wave.  These  two  cases  are  further  discussed  in 
Section  8  (although  the  simulations  in  Section  8  are  not  set  up  to  illustrate  the 
dynamics  near  these  nontrivial  equilibrium  points). 

In  a  linear  stability  analysis  (not  shown),  the  system  (8)  is  neutrally  stable 
when  linearized  around  the  trivial  equilibrium  point  of  0.  On  the  other  hand, 
when  linearized  around  the  two  examples  of  the  nontrivial  equilibrium  points, 
the  system  is  unstable.  It  would  be  interesting  to  further  explore  these  cases  and 
their  physical  interpretations  in  the  future. 


=  0,  a 

(27a) 

(27b) 

(27c) 


6.3  Derivation 

The  derivation  of  (24)  starts  by  assuming  a  leading  order  barotropic  component 
of  the  form 

Si  =  -4==0(ri,T2)ei*T  +  C.C.,  (29) 

V  -C/T 

where  0 t  =  k^x'  +  ,  and  is  the  barotropic  energy  unit 

Et  =  2  YL2, 

and  leading  order  baroclinic  variables  of  the  form 

m  =  a1ei6lf1  +  a2ei02r2  +  C.C., 


(30) 
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where  6j  =  kjx'  +uijt',  and  r  represents  a  right  eigenvector  of  the  linear  system 
with  components 


rj  -  (Ki,j ,  Rij ,  Qij ,  Aij ,  iC2)  ij ,  qt2')  1:j) ,  j  —  1,2. 


Here  the  arbitrary  complex  phase  factor  of  each  eigenvector  is  chosen  so  that 
K ij  is  a  real  number,  and  the  eigenvectors  fj  are  normalized  with  respect  to 
the  energy  as 

Ej=rfjH?j  =  1,  j  =  1,2,  (31) 

where  H  is  the  Hessian  matrix  of  the  conserved  energy  for  the  linear  baroclinic 
system. 

Next,  two  auxiliary  problems  must  be  solved,  (C.5)  and  (C.6)  of  Appendix  C, 
one  for  each  of  the  two  long  time  scales,  T\  =  St'  and  T%  =  52t' .  The  solution  of 
the  auxiliary  problems  is  the  key  step  in  the  method  of  multiscale  asymptotics 
in  order  to  suppress  secular  growth  and  guarantee  a  consistent  asymptotic  ex¬ 
pansion  of  the  variables  [13].  As  part  of  the  second  auxiliary  problem,  the  three 
waves  must  satisfy  the  following  resonance  conditions: 


fci  +  &2  +  =  0,  (32a) 

to  1  T  cd2  T  k-’T  =  0-  (32b) 


Further  details  of  derivation  are  presented  in  Appendix  D. 


7  Validation  study:  the  wavenumber-2  MJO 
mode 

In  this  section,  a  special  case  is  considered  in  order  to  explore  the  time  scales  of 
validity  of  the  ODE  system  in  (24).  If  only  a  single  baroclinic  mode  is  considered, 
without  a  barotropic  mode,  the  ODE  system  in  (24)  becomes  a  single  ODE: 

—  dT2<xi  +  id,4ai2ai*  +  id^ai  =  0.  (33) 

The  energy  conservation  of  this  single  mode  then  takes  the  form 

Jr(ai«i*)  =  0.  (34) 

This  special  case  can  then  be  compared  against  numerical  solutions  that  have 
been  presented  previously  elsewhere  [19],  which  will  be  regarded  as  the  ‘true’  so¬ 
lution  here.  The  specific  case  chosen  is  Case  U2  of  [19],  which  has  a  wavenumber- 
2  MJO  eigenmode  as  the  initial  condition  for  a  nonlinear  simulation. 
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A  few  technical  details  for  consistency  with  [19]  are  the  following.  The  coef¬ 
ficient  Ci  in  (A.l)  is  taken  to  be  1  for  consistency  with  the  equation  At  =  TQA 
used  in  [19],  Also,  since  the  model  of  [19]  does  not  include  the  dispersive  effects 
of  t/1)  (although  it  does  include  other  dispersive  effects),  the  coefficient  d§  in 
the  ODE  system  (33)  is  set  to  0.  As  in  [19],  here  the  computational  domain  is 
[0,  X],  where  X  =  40,  000  km  in  the  circumference  of  Earth  at  the  equator. 

For  the  validation  study,  three  different  values  of  5  are  chosen:  5  =  \/0.1, 
\/C)TT/2  and  \/0.1/4,  corresponding  to  three  different  RCE  states:  S2a  =  0.1, 
0.025  and  0.00625.  The  solutions  are  evolved  out  to  times  100,  200,  and  400 
days  respectively,  matching  the  long  wave  assumption.  Table  4  shows  the  rela¬ 
tive  error  between  the  asymptotic  solution  and  the  ‘true’  solution.  The  relative 
error  decreases  by  a  factor  of  2  as  5  is  decreased  by  a  factor  of  2,  indicating 
a  relative  error  of  0(5),  suggesting  an  absolute  error  of  0(53)  and  a  true  solu¬ 
tion  with  magnitude  of  0(52).  Figures  2  and  3  show  comparisons  between  the 
asymptotic  solutions  and  numerical  (‘true’)  solutions  for  5  =  x/OT  and  x/O.1/2 
(with  5  =  vOfL/4  not  shown  since  it  is  indistinguishable  from  the  ‘true’  solu¬ 
tion).  In  Figure  2,  the  agreement  is  good  out  to  100  days,  although  some  of  the 
peaks  in  convective  activity  HA  are  not  captured  by  the  asymptotic  solution. 
In  Figure  3,  the  agreement  is  excellent  out  to  200  days. 

As  a  further  validity  test  of  the  asymptotic  model,  the  same  comparisons 
have  been  repeated  (not  shown)  with  the  addition  of  the  second  order  cor¬ 
rections.  More  specifically,  whereas  the  asymptotic  solution  IFasym  in  Table  4 
was  defined  as  IFasym  =  52Wi ,  the  comparisons  have  been  repeated  using 
Wasym  =  52W\  +  53W2,  where  the  second  order  correction  W2  is  described 
in  more  detail  in  Appendix  C  and  D.  I11  these  new  tests,  the  values  of  Table  4 
change  to  0.2685,  0.0627,  0.0183.  These  values  show  a  decrease  in  relative  error 
of  a  factor  of  4  as  5  is  decreased  by  a  factor  of  2,  indicating  a  relative  error  of 
0(52),  in  line  with  an  absolute  error  of  0(<54)  and  a  true  solution  with  magnitude 
of  0(52). 

In  summary,  the  asymptotic  solutions  have  significant  accuracy  on  time 
scales  of  100  days  or  longer,  which  are  roughly  the  time  scales  for  application  to 
tropical-extratropical  interactions  and  MJO  initiation  and  termination. 
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52a  =0.1  at  100  day 

<5 2a  =0.025  at  200  day 

52a  =0.00625  at  400  day 

||WaSym-  Wn||/||Wn|| 

0.3890 

0.1657 

0.0770 

Table  4.  The  relative  difference  between  numerical  solution  Wn  (regarded  as  the  'true'  solu¬ 
tion  here)  and  the  asymptotic  solution  (lTasym)  whose  amplitude  cxj1(T2)  is  governed  by  the 
ODE  (33).  Three  different  values  of  S  are  considered  along  with  three  different  correspond¬ 
ing  times. 


d2 

dg 

C?4 

de 

dy 

ds 

dg 

MRB 

-0.45 

— 8.3e-3 

— 5.2e-2 

1.7e-3 

-0.28 

— 2.4e-2 

1.37 

0.29 

MKB 

0.45 

— 3.2e-2 

-0.11 

6.0e-3 

0.11 

7.8e-3 

— 8.5e-7 

—  7.6e-2 

Table  5.  Coefficients  dj  for  the  3-wave  interaction  ODEs  (24)  for  two  cases:  the  MJO,  dry 
baroclinic  Rossby,  and  barotropic  Rossby  wave  interactions  (MRB)  (top  row)  and  the  MJO, 
dry  Kelvin,  and  barotropic  Rossby  wave  interactions  (MKB)  (bottom  row). 


8  Three-wave  interactions 

Here  two  types  of  three-wave  interactions  are  considered  as  shown  in  figure  4: 

(1) :  MJO-dry  baroclinic  Rossby-barotropic  Rossby  waves  (MRB) 

^MJo(l)  +  Wr(- 3)  +  WT,yRi2,000  km(2)  =  0 

(35) 

(2) :  MJO-dry  Kelvin-barotropic  Rossby  waves  (MKB) 

wMJo(  —  3)  +  Wk(2)  +  ^T,y r*3,000  km(l)  =  0 

For  the  acronym  MRB,  the  M  stands  for  MJO,  the  R  stands  for  equatorial 
Rossby  wave,  and  the  B  stands  for  barotropic  Rossby  wave;  and  for  MKB, 
the  M  and  B  are  the  same,  and  the  K  stands  for  equatorial  Kelvin  wave.  The 
meridional  wavelength  2y  of  the  barotropic  Rossby  wave  is  chosen  so  that  the 
resonance  conditions  (32)  are  exactly  satisfied.  The  numerical  results  here  are 
focused  on  the  ODE  system  (24),  with  wave  amplitudes  (3,  a±,  and  «2,  where 
the  energy  exchange  between  three  waves  is  of  the  most  interest.  The  coefficients 
in  (24)  for  the  two  cases  are  in  table  5.  The  RCE  a,  is  chosen  to  be  a  =  0.1  so 
that  the  small  parameter  5  =  x/cTT. 
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8.1  MJO-dry  baroclinic  Rossby-barotropic  Rossby  waves 
(MRB) 

In  this  case,  a±  and  <*2  are  the  amplitudes  for  the  MJO  and  dry  baroclinic  Rossby 
waves,  respectively.  To  consider  MJO  initiation,  a\  is  set  to  be  zero  for  the  initial 
condition  in  the  ODE,  and  the  other  waves  have  amplitudes  of  |a2(t=o)l  = 
|/3(t= o)|  =  1.  The  results  of  the  numerical  simulation  are  shown  in  Figure  5. 
The  waves  exchange  energy  periodically,  with  period  of  about  60  days,  roughly 
the  time  scale  for  MJO  initiation  and  termination  in  nature.  At  early  times, 
the  MJO  amplitude  grows  (i.e.,  the  MJO  initiates)  by  extracting  energy  from 
the  other  two  waves,  mainly  from  the  dry  baroclinic  Rossby  wave.  Termination 
of  the  MJO  then  follows  from  a  transfer  of  energy  back  to  the  dry  baroclinic 
Rossby  wave. 

To  explore  the  sensitivity  of  these  results,  different  initial  conditions  are  cho¬ 
sen,  and  the  results  are  in  figure  6.  The  initial  amplitudes  for  baroclinic  Rossby 
wave  and  barotropic  Rossby  waves  are  (a)  |c*2(t=o)l  =  0-5,|/9(t=o) I  =  1;  (b) 
|<*2(t=o)l  =  1)  |/?(t=o)l  =  0.5;  (c)  |a2(t=o)i  =  0.25,  |/3(t=0)|  =  1;  (d)  |a2(t=o)l  =  1; 
|/3(t=0)l  =  0.25.  The  simulations  are  performed  with  different  initial  phases  of 
«2  and  /3,  but  this  has  no  effect  on  the  evolution  of  the  amplitudes.  In  the  sim¬ 
ulations  in  Figure  6,  it  can  be  seen  that  the  amount  of  energy  extracted  by  the 
MJO  varies  roughly  proportionally  to  the  amount  of  energy  present  initially  in 
the  dry  baroclinic  Rossby  wave  and  barotropic  Rossby  wave.  On  the  other  hand, 
the  time  scale  for  the  energy  exchange  changes  very  little  among  these  cases.  In 
all  cases,  the  energy  of  the  barotropic  Rossby  wave  stays  essentially  constant  and 
is  not  transferred  to  the  MJO.  These  results  indicate  that  the  eventual  strength 
of  the  MJO  depends  on  the  amplitude  of  the  barotropic  Rossby  waves,  but 
the  MJO  acquires  its  energy  from  the  baroclinic  rather  than  barotropic  Rossby 
waves. 


8.2  MJO-dry  Kelvin-barotropic  Rossby  waves  (MKB) 

In  this  case,  a\  and  «2  are  the  amplitudes  for  the  MJO  and  dry  baroclinic  Kelvin 
waves,  respectively.  Again  considering  MJO  initiation,  the  initial  conditions  are 
{*1=0  and  |o!2(t=o)l  =  |/3(t=o) I  =  1,  and  the  results  are  shown  in  Figure  7. 
While  the  solution  is  similar  in  character  to  the  MRB  case,  two  main  differences 
can  be  seen.  First,  the  temporal  period  of  roughly  120  days  is  twice  that  of 
the  MRB  case  but  still  within  the  range  of  time  scales  for  MJO  initiation  and 
termination  in  nature.  Secondly,  both  the  dry  Kelvin  wave  and  the  barotropic 
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Rossby  waves  contribute  to  the  growth  of  MJO,  although  the  energy  exchange 
with  the  barotropic  Rossby  wave  is  somewhat  small. 

Finally,  the  sensitivity  of  the  results  is  illustrated  in  Figure  8,  which  shows 
numerical  results  with  different  initial  amplitudes.  As  in  the  MRB  case,  the 
amount  of  energy  extracted  by  the  MJO  varies  roughly  proportionally  to  the 
amount  of  energy  present  initially  in  the  other  two  modes. 


9  Summary  and  conclusions 

A  new  model  has  been  proposed  here  to  describe  MJO  initiation  and  termina¬ 
tion  and  tropical-extratropical  interactions.  The  model  involves  the  integration 
of  the  barotropic  and  equatorial  baroclinic  modes  together  with  moisture  and 
convective  activity  envelope,  all  interacting  together  with  a  conserved  energy. 
Using  the  method  of  multiscale  asymptotics  with  multiple  time  scales,  simpli¬ 
fied  asymptotic  equations  were  derived  for  the  resonant  interaction  of  tropi¬ 
cal  and  extratropical  waves.  The  simplified  model  is  an  ODE  system  for  wave 
amplitudes,  including  quadratic  nonlinearity  as  in  many  traditional  three-wave- 
resonance  equations,  and  also  including  cubic  self-interaction  terms  that  arise 
from  nonlinear  interactions  of  water  vapor  and  convective  activity. 

Example  simulations  of  the  ODE  system  are  shown  to  illustrate  some  cases 
of  MJO  initiation.  In  this  model,  the  MJO  is  shown  to  initiate  mainly  by  extract¬ 
ing  energy  from  either  the  dry  Kelvin  wave  or  the  dry  equatorial  baroclinic  Rossy 
wave.  While  the  MJO  extracts  a  smaller  amount  of  energy  from  the  barotropic 
Rossby  waves,  the  barotropic  Rossby  waves  are  essential  for  MJO  initiation, 
and  the  strength  of  the  ensuing  MJO  depends  on  the  strength  of  the  barotropic 
Rossby  waves  present.  In  this  way,  the  barotropic  Rossby  wave  acts  as  a  catalyst 
for  the  interaction  of  the  other  two  waves.  In  other  scenarios,  it  is  possible  that 
the  barotropic  Rossby  wave  may  lose  or  gain  energy  through  wave  interactions 
if  more  realistic  setups  are  used.  For  example,  (i)  zonally  varying  climatological 
mean  state,  and  (ii)  shear  could  affect  the  interaction  mechanisms.  In  the  fu¬ 
ture,  it  would  also  be  interesting  to  include  the  effects  of  additional  meridional 
modes  that  are  asymmetric  with  respect  to  the  equator  in  order  to  model  the 
Boreal  summer  MJO  and  monsoon  intraseasonal  variability  [35].  The  authors 
are  currently  pursuing  these  issues  and  will  report  on  them  in  the  near  future. 

The  simplified  models  suggest  specific  classes  of  waves  that  could  be  inves¬ 
tigated  in  observations  among  the  plethora  of  possible  waves  and  wave  inter¬ 
action  scenarios.  More  specifically,  one  can  imagine  wave  interactions  among 
many  different  types  of  waves.  Here,  in  this  paper,  we  have  identified  specific 
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classes  of  three-wave  interactions  involving  the  MJO.  For  example,  one  class 
involved  a  wavenumber-1  MJO,  wavenumber-3  dry  baroclinic  Rossby  wave,  and 
wavenumber-2  barotropic  Rossby  wave.  It  would  be  interesting  to  identify  this 
specific  class  of  waves  in  observational  data  and  investigate  their  interactions. 
Such  an  investigation  and  comparison  with  the  model  could  be  challenging  for 
many  reasons;  for  example,  as  described  above,  the  mean  state  of  the  model 
version  here  does  not  include  a  zonally  varying  climatological  mean  state  nor 
wind  shear,  which  may  influence  the  wave  interactions. 


A  The  explicit  formulation  of  the  meridional 
truncated  system 

Applying  the  meridional  truncation  of  (14)  to  the  long- wave-scaled  system  of 
(6)  leads  to  the  truncated  system: 


L2YBt, 

-  YBX,  =Ct  +  52  (YBx>x>t,  +  VT) 

(A.la) 

- 

l(0)x,  +  „(!)  +  l=5Ha{  0)  =  Cm 

v  2 

(A.lb) 

l{2\  ■ 

-  —  Cn 2) 

(A.lc) 

r(°)t, 

+ 

r(0)*'  +  =  Cr(  o) 

(A.ld) 

r^t, 

+ 

r^2\’  -  V'2'd(1)  =  Cr( 2) 

(A.le) 

s/2rW  . 

-i(°)  =  52(-vW'tl+Vvw) 

(A.lf) 

+ 

^(r(0)^  -  1{0K')  +  ^U(1)'  +  SHa^  =  Cq(0) 

(A.lg) 

?(2V 

+ 

^(r(2)a;'  -  l(2)x')  ~  Qv^'  =  Cq( 2) 

(A.lh) 

a<°>. 

— 

rag(0>  =  cira(oyo)  +  Ci0ra(°)g(2) 

(A.li) 
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which  is  written  in  vector  form  in  (19).  Here  C  and  V  are  bilinear  terms  at 
different  orders: 


CT  =  -|(r( 0)  _  l( 0))2,  _  |  J(r(0)  _  i(0))(r(2)  _  Z(2)} 

_  £i(r(2)  _  /(2))2  _  (rC0)  _  ;(0))t,(l)'  _  (r(2)  _  j(2))w(l)' 

2  ,x  2^2  2V2 


(A. 2a) 


C,(«  =  ~c2BxZ  -  c&BxdW  -  2 c2bZx,  -  c3B^2>xr 

+  c2(r(0)  -  Z^)^,  +  — (r(2)  -  Z^2))Bx/  -  (A.2b) 

2  v2 

C;(2,  =  - c9Bx,l (°)  -  c4Sx^(2)  -  C3BZ(0)X'  -  2 c4BZ(2)^ 

+  -(r®  -  Z(0))HX/  +  c4(r(2)  -  Z(2))HX/  -  (A.2c) 

2  V2 

Cr(0)  =  -c2Bxirw  -  c8Bxrr(2)  -  2 c2Sr(0)x<  -  c3Sr(2)x/ 

-  c2(r(°)  -  j(°))Bxr  -  ^(r<2>  -  ?(2))Bx,  +  v^' B  (A.2d) 

Cr(2)  =  -c9Hx/r(0)  -  c4Sx/r(2)  -  c3Br(0)l(  -  2 c4Br(2)^ 

-  7r(r(0)  “  ^(0))^'  -  c4(^(2)  -  1{2))BX'  +  -^=vw'b  (A.2e) 

2  v  2 

C,(0,  =  -c2Bx,qW  -  2 c2Bq^x,  -  c&Bx,qW  -  c3Bq^x,  (A.2f) 

Cq( 2)  =  -c4Bx/g(2)  -  2 c4Bg(2)x,  -  c9Sx/g(0)  -  c3Bg(0)x,  (A.2g) 


and 


= 


' 


\[2L2 


(r(0)  _  Z(0))V(D' 


C6 


V2L2 


(r(2)_z(2))v(l)'' 


+  20^ 

(A. 3a) 


2?„(i)  =  -c7.Bay^!(1),  -  2c75u(1)x,  +  2c7t)(1),Bx 


+  2=  [gr(r(0)  -  l{0))BX'x'  +  §(^(2)  -  ^(2))^'*.']  (A. 3b) 

where  the  coefficients  CjS  are  defined  by  parabolic  cylinder  functions: 


Cl=fZ*jdy, 

c3  =  —L  fZ  $n$2cos(Ly)dy, 
C5  =  -L2  fZ  sin(Ly)dy, 
c7  =  -4fZo  cos(Ly)dy, 
c9  =  I-00  $Q$2  sm(Ly)dy, 


C2  =  ~k  LZ  Hcos(Ly)dy , 

C4  =  -f  fZ  $2  cos(Ly)dy, 
c6  =  -L2  fZ  $i$2  sin(Ly)dy, 

08  =  f-o o  $0$,2  sin(L^)dy> 

cio  =  fZ  $o$2d y, 


(A.4) 


Notice  that  the  equality  cs  +  c9  =  c3  helps  to  conduct  the  energy  conservation. 
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The  truncated  system  in  (A.l)  conserves  a  total  energy  if  the  term  a(°)qi(2) 
is  neglected  in  (A.li): 

x 

£=\f  2YL2(B2x  +  B2)  +  |(°)2  +  l^2  +  r(°)2  +  r^2  +  -  alog(a<°>  +  a)) 

o 


1 

+  <0(1  -  Q) 


Q_ 

72 


2 


(r(°)  +Z(°)) 


+ 


72 


2 


(r(2)  +Z(2)) 


da: 

(A.5) 


For  this  reason,  the  term  £j(°)g(2)  in  (A.li)  will  be  neglected  throughout  the 
present  paper. 


B  Asymptotic  expansion  of  the  truncated 
system 


Here  an  explicit  expansion  is  presented  of  the  truncated  system  from  (A.l), 
which  was  written  in  vector  form  in  (19).  Using  an  asymptotic  expansion  to  the 
third  order,  the  ansatz  is: 

B  =  52B1  +  53B2  +  54B3  +  0(7) 

(/w  ^V((V(2\71)/,70),72)) 

=  7  (Z(0)  r ,  Z(2) ! ,  r (0) ! ,  r (2) ! ,  w (1)  r ,  q(0)  x ,  q(2\) 

+  63(l(0)2,l(2)2,  r(°)2,  r(2)2,v(1)2,  q(0)2,  q{2\) 

+  7(/W3,  r<°>3,  r®3,  «W3)  70)3, 72)3)  +  0(7)  (B.la) 

a(°)  =  5a(0)1  +  ^a(0)2  +  53a(0)3  +  q(5  4)  (B.lb) 

All  variables  are  assumed  to  depend  on  three  time  scales:  t'  and  T\  =  St1  and 
T2  =  52t'. 

To  simplify  notation,  define  Un  =  (Z(0)„,  l^n,r^n,  r^n,  v(1)n,  q(0)n,  q(2)n,  a(0)n), 
and  let  £{y  denote  the  spatial  linear  operator. 

Using  this  ansatz,  the  truncated  system  (A.l),  or  its  vector  formulation 
(19),  can  be  asymptotically  expanded  over  three  orders  of  magnitude.  The  first 
order  system  is 


L2YBw  -  YBlx,  =  0, 
A fUlt'  +  £vUi  =  0, 


(B.2a) 

(B.2b) 
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The  second  order  system  is 


L2YB2t>  -  YB2x,  =  -L2YB1Tl, 

MUW  +  C-ijU 2  =  —U\Tt  +  Fq  2, 

The  third  order  system  is 

L2YB3t,  -  YB:ix,  =  -L2YB1T2  -  L2YB2Tl+YBlx,x,t,  +C1T, 


(B.3a) 

(B.3b) 

(B.4a) 

(B.4b) 


MVw  +  3  =  ~U\t2  ~  +  Fjj3. 

Here  J\T  =  diag(l,  1, 1, 1, 0, 1, 1, 1)  is  the  8x8  matrix  just  to  eliminate  dt>  for  the 
variable.  Forcing  terms  Fjj2  and  Fq3  are  quantities  from  the  lower  order: 

Fjj  2  =  (0,0,0,0,0,0,0,cira(°)1gW1),  (B.5a) 

Fjj3  =  (Cu(o)  ,Cu(2), 

Cir(o) ,  Cir(2) ,  — i/1)  it' ,  Clq( o) ,  Clg(2) ,  ciT(a(0)  \q'a\  +  i))- 

(B.5b) 

Here  C\  denotes  the  bilinear  operations  are  performed  on  U\.  Like  the  2D  asymp¬ 
totic  expansion,  system  (B.2)-(B.4)  carries  similar  approximated  energy  conser¬ 
vation  that 

=  o(54).  (B.6) 


C  Details  of  the  auxiliary  problem 

In  this  appendix,  the  auxiliary  problem  is  written  explicitly  in  terms  of  the  new 
variables  W  =  ( K ,  R,  Q ,  A,  v^\  x),  where  the  change  of  variables  was  described 
abstractly  in  section  4. 


C.l  Reformulation  of  the  auxiliary  problem 

Equation  (21)  can  be  concretely  written  as 

Kf  +  Kx  4 — —HA  =  Fr( o) 

v  2 

Rt  +  V2\x  ~  \f2v^  +  HA  =  '/2Fl(0)  +  2Fr(2) 

1  ^  Q 

2Xx)  +  V2 

,  At  -YaQ  =  Fa( o) 


Qt  +  %  {Kx  2V2 


(C.la) 

(C.lb) 


Rx  +  oXx  )  +  +  HA  —  F  t o)  (C.lc) 


(C.ld) 
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0-v(1)t  +  x  =  ^u(i 


) 

Xt  +  —  3'1'*'1'  —  =  —  -^(O)  + 


(j(2) 


(C.le) 

(C.lf) 


Z(2)t  -  1{2K  =  Fti 2) 

?(2)* +  +  xXx' _  z(2)*')  _  =  F<- 

Tlie  equations  (C.le)-(C.lf)  are  used  to  eliminate  V1)  and  x  in  the  other  equa¬ 


(C-lg) 

(C.lh) 


tions: 


X  =  Kd)  (C.2a) 

vW  =  (**  -  HA)  +  i  (Fo(1)it  +  Fm  -  V2 FrW)  (C.2b) 

so  that  the  system  becomes 


Kt  4-  Kx  4 — —HA  —  Fr(o) 
v  2 


1 


4  - 


Ri  +  —  Ra,  ~b  —  H  A  — 


3  3 

IV 

3 


4\/2  4  r  \/2 

—^—Fn o)  +  —Fr(2)  -  V2F^(i)>x  +  —Fv( i)>t 


(C.3a) 


(C.3b) 


+  ~^Fr<-2')  +  Fq(o)  -  -^j=FvW  x  -  ~^J^FV( i)>t(C.3c) 


,  —  TaQ  —  -Fa(o) 


(C.3d) 


l(2)t-l(2K-  =  F, 


q{2\- 


-^rx-^2K  +  ^a  = 
12V2  V2  3V2 

P  1  *3  p  V2Q  Q 

t,(2)  +  y-rj(o) - 7^—rA2)  - 


(C.3e) 


Q 


~  (C.3f) 


When  jV  =  0,  the  system  (C.3),  (C.3a)-(C.3d)  forms  the  KRQA  system  as 
in  [18],  while  the  equation  (C.3e)  for  Z^2)  is  independent,  and  equation  (C.3f)  for 
is  slaved  by  KRQA  system  and  Z^2). 


C.2  Forcing  vector  for  the  auxiliary  problem 

Recall  the  earlier  presentation  of  auxiliary  problems  in  (B.2)-(B.4)  in  terms 
of  the  variables  U.  In  the  transformation  from  U  to  W  that  was  described 
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abstractly  in  section  4,  the  forcing  vector  Fq  is  transformed  into  forcing  vector 
Fyfr.  The  transformation  of  the  forcing  vector  is  now  presented  explicitly  for 
each  of  the  cases  in  (B.2)-(B.4).  Using  the  notation 

B  =  52Bi  +  53B2  +  54B3  +  0(S5), 


and 

W  =  52Wi  +  53W2  +  54W3  +  0(<55), 

the  leading  order  system  is 

L2YBlyyt,  -  YBlxr  =  0  (C.4a) 

Ww  +  C^W i  =  0  (C.4b) 

The  second  order  system  is 

L2YB2yyt'  —  Y  B2xi  =  —Biyyx1,  (C.5a) 

w2t>  +  C^W2  =  -W1Tl  +  Frf2,  (C.5b) 

and  the  third  order  system  is 

L2Y  B3yyt'  —  Y  B3xi  =  —BiyyT2  —  B2yyTl  +  Bixixit/  +  Cit  (C.6a) 
w3t>  +  C,^W3  =  -W1T2  -  W2Tl  +  U#3.  (C.6b) 

Here 

J?#2  =  (O.O.O.ciTAiQi,  0,0)  (C.7a) 

F\y3  =  (^lr(°)>  — ^ ,C0)  3^lU2>  +  V/2«^1')l x't'  ~ 

Q  n  i  Q r  i  r  Q  (l)  ,  Q  „.(i) 

“37iCli(0>  +  JClr(2)+Cl^0)  +  ^v  lx't,  +  ^v  14,4,1 
cir(AiQ2  +  A2Qi)  +Cla(  o),CU(2), 

Clq  (2)  +  -jj-CllW - -^~Clr(  2)  +  (C.7b) 

Note  that  for  the  first  and  second  order,  xi  =  X'2  =  0,  and 
v{1)n=  ^-j=(Rnx' -  HAn)  ,  n  =  1,2. 


(C.8) 
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D  Details  for  multiscale  analysis 

D.l  Solving  linear  system  with  forcing 

Here  a  brief  summary  for  solving  linear  system  with  forcing  is  provided,  which  is 
a  guideline  of  the  multiscale  analysis.  Readers  may  find  information  from  other 
references,  e.g.,  [13].  Let  W  solve  the  linear  system  with  forcing: 

dt-W+  C^W  =  Pe^kx-ct'\w\t'=o  =  0.  (D.l) 

Next,  assume  that  W  is  the  superposition  of  eigenmodes: 

=  J2as(t')rs(k)ei(kx'~u;i>{k)t')  ■  (D.2) 

S 

Two  cases  might  happen: 


1.  If  us(k)  ±  c,  then  as  = 

2.  If  uis(k)  =  c,  then  as  =  t  ( ls  ■  Fj. 

At  the  second  case,  when  !sT/0,  linear  growth  in  time  will  happen.  Therefore 
the  solution  W  is  bounded  if  either  1)  ujs(k)  ^  s  for  all  s;  or  2)  ws(fc)  =  c  and 
ls  ■  F  =  0.  The  above  is  the  principal  of  the  multi-scale  analysis. 

D.2  Solutions  for  the  second  order  equations 

At  the  second  order,  the  only  nonlinear  term  is  in  the  A  equation  in  (C.5).  The 
barotropic  and  baroclinic  modes  are  still  decoupled.  They  are  studied  separately. 
For  the  baroclinic  wave,  the  only  nonlinear  term  does  not  generate  resonance 
condition  between  different  eigenmodes.  Therefore,  considering  one  mode  is  suf¬ 
ficient  to  study  solutions  for  the  second  order  equations. 

Assume  that  the  leading  order  solution  has  a  form  of 

Wi  =  a1(T1,T2)ei01fjl  +  C.C. 

Here  9 i  =  k\x '  —  wi(fe1)t/  is  the  phase,  and 

fj'i  =  (ATi,i,  -Ri,i,  Qi,i,  ^i,i,  9^1,1.) 
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is  the  eigenvector  for  the  k  =  k\  eigenmode  of  certain  wave  type  (e.g.,  Kelvin, 
MJO,  moist  Rossby,  dry  Rossby). 

Also  assume  that  the  leading  order  baroclinic  wind  is 

Wi  =  a\ei6lri  +  a2ei62r2  +  C.C., 


where  9j  =  kjX1  +  uijt ' ,  and 

?j  =  j  =  1,2. 

Here  the  eigenvectors  are  chosen  so  that  Kijs  are  real  numbers.  The  eigenvector 
fj  are  normalized  by  their  energy  unit,  so  that 

Ej=?)nr0  =  1,  j  =  1,2,  (D.3) 

where  H  is  the  Hessian  matrix  of  the  conserved  energy  for  the  linear  system. 

At  the  second  order,  the  forcing  terms  come  from  the  leading  order  (C.5). 
The  solution  consists  of  two  parts:  the  homogeneous  solution  and  the  nonho- 
mogeneous  solution  from  the  forcing  terms.  The  homogeneous  solution  are  just 
linear  eigenmodes  which  is  considered  to  be  absorbed  into  the  leading  order. 
Therefore  the  second  order  solution  is  completely  determined  by  forcing  terms: 


F  =  —a\dTlWi  -  a\dTlW*  + 


0 

0 

0 

^3 

0 

0 


(D.4) 


where 

F3  =  ci T  (^ai2ei2Sl  Aij1Qij1  +  a^A^Q^)  +  C.C. 
The  right-hand-side  of  the  second  order  has  three  phases: 

e±ie\e±i29\  and  e°. 


1.  Phase  e±lBl 

These  two  phases  generate  secular  growth.  The  phase  e*®1  requires  that 

^i(fci)  '  ^(fci)  =  d^ai  =  0. 

The  same  thing  applies  for  phase  e-*®1,  so  that  d^oti  =  0. 
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2.  Phase  e±l2e° 

This  part  comes  from  the  Q-A  nonlinear  term.  Because  it  does  not  generate 
secular  growth,  the  system  can  be  solved  and  the  solution  is 


W, 


2(2fci  ,t',T2) 


n{ 2fcl)(e-2i“1('=DtI  -  e-v^t)(V{2kl)  -  2 iL01{kl)iy 


•C(2fei)-P1(2fei) 


(D.5) 


3.  Phase  e° 

This  part  also  stems  from  the  Q-A  nonlinear  term.  Although  it  has  0- 
frequencies,  which  may  resonant  with  the  MJO,  moist  Rossby  mode,  0- 
frequency  mode,  fast  west-propagating  l ^  mode,  but  it  does  not  happen 
because 


^MJO(O)  ‘  *(0)  -  ^mR(O)  '  F{0)  ~  ^Ofr(O)  ‘  F{0)  ~  h<-2) ( 0)  '  F(0)  ~ 

which  is  also  stated  in  [13].  The  phase  e°  part  have  the  contribution  to  the 
solution  as 


( 


W, 


2(fe=0,  t',T2) 


=  n 


(0) 


diag 


V 


1 

ilVK(0) 

0 

0 

0 

1 

“R(O) 

0 


\ 


(x_e-®(o,t')£{o)J;(o) 


/ 


One  remark  is  that  because  the  leading  order  as  forcing  terms  does  not  have  T\ 
dependence,  the  second  order  solution,  is  also  independent  of  Tf,  i.e., 


d 

dT\ 


#2  =  0. 


D.3  Resonance  condition  for  the  third  order  equation 

At  the  third  order,  there  are  nonlinear  terms  mixing  baroclinic  and  barotropic 
terms,  together  with  the  nonlinear  q-a  interaction. 
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To  simplify  mathematical  calculations,  the  W\  variables  are  transformed 
back  to  the  U\  variables,  where  the  Fourier  coefficients  can  be  written  as 

=  ^yl^tdr  r(2)i,i  =  4^1,  h 
=  Qi,j>  a^i,j  = 

i’^ij  =  3^2  t-^tj  —  HA\ j)  ,  j  =  1,2. 

Here  the  three  groups  of  terms  in  (24)  are  explained  separately: 

1.  Cubic  terms. 

The  cubic  terms  in  (24)  are  associated  with  the  product  of  seconder  order 
and  the  first  order  from  the  nonlinear  q-a  interaction.  These  terms  are  self¬ 
interaction  -  they  do  not  interact  with  other  eigenmodes. 

Let  (A'2,  R-2,  Q2,  A2, 1^2,  9^2)  be  the  solution  obtained  in  Section  D.2. 
The  term  ciT(AiQ2  +  A2Q1)  contributes  to  the  resonant  condition  by 
the  nonlinear  interactions  between  the  first  and  second  order.  The  lead¬ 
ing  order  has  phases  e±l(klX  )j  whereas  the  second  order  has  phases 

e±i(2fci®'  — 2wi(fci)t')  e±i(i2fcix'— w(2fci)f')  e±iuiK( 0)t'  e±iuiR(0)t'  an(p  g0 

product  of  the  phases  -2wi(fci)t )  an{j  eo  jn  the  second  order  to¬ 
gether  with  the  phases  e±l(fclX  )  at  the  leading  order  will  contribute 

to  secular  growth. 

The  coefficient  idi  can  be  written  as 

d4  =  ft  •  (0,0,0,  d44, 0,0)  (D.7) 

where 


CL14 


—  ^l,l(-fci)Q2,l(2fci)  +  Ql,l(-fei)^2,l(2fci) 
+  J4l,l(fci)(32,l(0)  +  Ol,l(fei)^2,l(0) 


(D.8) 


where  h  is  the  left  eigenvector  of  the  choice  of  eigenmode  at  the  leading 
order,  and  wWpi  is  in  terms  of  1  and  Aqi  from  (C.8),  that 


And  idf  in  the  equation  for  02,  same  calculations  are  applied,  but  starting 
from  another  eigenmode. 

2.  Linear  terms 

Linear  terms  come  from  dispersive  effect  in  long- wave  scaled  system.  These 
are  also  self-interaction  terms. 
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For  the  barotropic  part,  the  coefficient  for  the  linear  term  *g?2  writes 

*^2  =  T9*^T2WT- 

L- 


For  the  baroclinic  part,  the  coefficient  ids  writes 

ids  =  h  •  (0,  ds2,  ds3,0, 0,  dse), 


where 


^  a. 

ds2  =  H — — uq2^1)  14, 

ds 3  = 

dsG  = 


(D.9) 

(D.10) 

(D.ll) 

(D.12) 


3.  Quadratic  terms 

The  quadratic  terms  are  associated  with  nonlinear  baroclinic-barotropic 
interactions. 

In  (24),  coefficient  ids  can  be  written  as 


ids  =  ~ 


1 


VEt^1 


(D-13) 


where 


d3i  =  -  i{k\  +  fc2)2c2(r(°)i)i  -  Z(°)iii)(r(°)ii2  -  K°) ij2) 

-  i(ki  +  fc2)2c4(r(2)i)i  -  l^\,i){A2\}2  -  Z(2)i,2) 

-  ic3(fci  +  fc2)(r(A°)iii  -  Z(°)i,i)(r(2)ij2  -  l&) i;2) 

-  ic3(k i  +  fc2)(r(°l  1,2  -  l(°)ii2)(r(2)i)i  -  Z(2)i,i) 

(r(°)i  i  _  |(o)lil)ti('1ll!2  +  (r(°)i)2  -  Z^i, 2)1^1, l 

"(r(”2)M  _  i(2)ljl)t,(i)  2  +  (r(2)12  _ 


_C5 

V2  L 

_C6 

L 


(D.14) 


and  for  baroclinic  mode,  coefficient  ids  can  be  written  as 


~id6  '  (Km 2^~^~Km2  +  ^Km2,~^j^Flm2  +  §Kma  +  Km 


9'UV 


Q  . 

Ft(2)„,  —  Ftm2  ~ 


a(°) 2  ’  ^ l(2^2  ’ 


r(2). 


+  ^9(2)2) 


3 


(D.15) 
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where  l±  is  the  left  eigenvector  for  the  j  =  1  baroclinic  mode,  and  F 
are  Fourier  coefficients  of  interactions  between  j  =  2  baroclinic  mode  and 
barotropic  wind,  They  corresponds  to  bilinear  terms  Cs  in  (A. 2).  One  ex¬ 
ample  is  given  for  F^o)'. 


-F/( 0)2  =  —  ifcx  ^C2/^l,2  +  Cal&  1,2^  —1^2  ^2c2i^l,2  +  C3ZW  1,2 

+  jfex  C2(r(°\t2  —  1^1,2)  +  —  1^1,2) - ^=V^l,2 

2  J  y/2 

(D.16) 


The  coefficient  dg  is  analogous  to  dg  by  swapping  j  =  1  and  j  =  2. 
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Fig.  2.  Comparisons  between  asymptotic  solutions  (left)  and  numerical  (‘true’)  solutions 
(right)  for  A  =  0.1  for  the  case  of  a  wavenumber-2  MJO  mode.  The  convective  activity  HA 
(top)  and  the  water  vapor  Q  (bottom)  are  shown  as  functions  of  x  and  t. 
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Fig.  3.  Comparisons  between  asymptotic  solutions  (left)  and  numerical  (‘true’)  solutions 
(right)  for  A  =  0.025  for  the  case  of  a  wavenumber-2  MJO  mode.  The  convective  activity 
HA  (top)  and  the  water  vapor  Q  (bottom)  are  shown  as  functions  of  x  and  t. 
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Fig.  4.  Resonance  condition  for  three-wave  interactions,  where  circles  indicate  the  particular 
wave  numbers  and  frequencies  that  leads  to  resonances.  Left:  MJO,  Kelvin  and  barotropic 
Rossby  wave;  right:  MJO,  Rossby  and  barotropic  Rossby  wave. 
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Fig.  5.  MJO— Rossby— barotropic  Rossby  wave  interactions.  Initial  condition:  ai  =  0, 

l“2(t=0)l  =  l/3(t=0)l  =  1. 
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Fig.  6.  MJO— Rossby— barotropic  Rossby  wave  interactions.  Initial  condition:  Qi  =  0.  (a) 

l«2(t=0)I  =  0.5,|/J(t=0)|  =  1;  (b)  |a2(t=0)l  =  1.  |j9(t=o)|  =  0.5;  (c)  |a2(t=0)|  =  0.25, 
l/3(t=o)l  =  1;  (d)  |«2(t=0)l  =  1.  l/3(t=o)l  =  0.25. 
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Fig.  7.  MJO— Kelvin— barotropic  Rossby  wave  interactions.  Initial  condition:  cri  =  0, 

l“2(t=0)l  =  l/3(t=0)l  =  !■ 
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Fig.  8.  MJO— Kelvin— barotropic  Rossby  wave  interactions.  Initial  condition:  qi  =  0.  (a) 

|«2(t=0)l  =  0.5,|/3(t=0) |  =  1;  (b)  |a2(t=o)|  =  1.  |/3(t=0)l  =  0.5;  (c)  |a2(t=0)|  =  0.25. 
|/3(t=0)|  =  1;  (d)  |a2(t=0)l  =  1,  l/3(t=0)l  =  0.25. 


